Prozen Disorder in a Driven System 
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We investigate the effects of quenched disorder on the universal properties of a randomly driven 
Ising lattice gas. The Hamiltonian fixed point of the pure system becomes unstable in the presence 
of a quenched local bias, giving rise to a new fixed point which controls a novel universality class. 
We determine the associated scaling forms of correlation and response functions, quoting critical 
exponents to two-loop order in an expansion around the upper critical dimension d c = 5. 
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The statistical mechanics of systems which settle into 
non-equilibrium steady states (NESS) has attracted con- 
siderable interest recently Q . Much attention focused on 
Ising-like systems subject to a driving force. These mod- 
els exhibit many intriguing features, such as anomalous 
correlations at all temperatures or order-disorder transi- 
tions belonging into various non-equilibrium universality 
classes [GJ . Yet their basic building blocks are designed to 
be as simple as possible, so that they are easily simulated 
and accessible by analytic methods. Thus, key character- 
istics of generic non-equilibrium behavior can be traced 
to specific model ingredients, a procedure which would 
be excessively difficult if one attempted to capture the 
full complexity of a real system. 

Microscopie models of driven systems typically involve 
nearest-neighbor particle-hole exchanges on a fully pe- 
riodic regular lattice, according to a sequential dynam- 
ics which is local and translation-invariant, in space and 
time. The hopping rates are controlled solely by (i) dif- 
ferences in internal energy, determined by a local Hamil- 
tonian, (ii) the coupling to a heat bath which enters 
through its temperature T, and (iii) a driving force. 
The latter acts as a dynamic perturbation which drives 
the system away from thermal equilibrium, into a non- 
Hamiltonian steady state which depends on the details 
of the dynamics. The prototype is a uniformly driven 
system B , motivated by the physics of fast ionic con- 
ductors B: Here, the usual Ising energetics is modi- 
fied by a uniform bias which favors (suppresses) jumps 
along (against) a specific lattice direction, thus inducing a 
global current through the system. Recast as a Langevin 
equation in continuous space and time, the theory can 
be analyzed using mean-field or renormalization group 
methods. Since, in principle, this procedure relies on 
coarse-graining the microscopie dynamics, the resulting 
equation is identical for all choices of local microscopie 
rules which respect the same set of global symmetries 
and conservation laws. In contrast, different universal 
behavior, controlled by a modified Langevin equation, is 
expected if the dynamics is nonlocal, in space or time. 



In this letter, we consider the effect of nonlocality in 
time. In any experimental realization of driven lattice 
gases, the presence of disorder, in the form of impurities 
or lattice defects, is to be expected. If such defects are 
effectively frozen on the time scales that govern the the 
order parameter dynamics, they introducé long (techni- 
cally infinite) correlation times, or memory effects, into 
the dynamics, which may lead to striking consequences 
for static and dynamic critical behavior. In equilibrium 
systems, the study of quenched disorder has a long his- 
tory pj], and some general features, such as the Harris 
criterion Q for the stability of pure fixed points, can be 
identified. For driven systems, given the absence of a free 
energy, it is not known whether any such simple criteria 
exist. We therefore present a case study of a specific 
driven system, as a step towards a more general explo- 
ration of the effects of quenched randomness on NESS. 

The choice of our model is motivated by two earlier 
studies on frozen disorder in the uniformly driven lattice 
gas, one simulational || and the other analytic 0, con- 
sidering, respectively, a small concentration of immobile 
impurities versus quenched randomness in the nearest- 
neighbor hopping rates |j. In both cases, thermody- 
namic observables are averages over different disorder re- 
alizations. However, the presence of a global current, in- 
duced by the uniform drive, leads to a fundamental diffi- 
culty: by necessity, the simulations are performed on a fi- 
nite periodic lattice, and before the system reaches steady 
state, individual particles have been driven through the 
lattice several times, encountering the same disorder con- 
figuration on each visit. Thus, the usual time averages 
include uncontrolled correlations in the randomness, in 
stark contrast to the physics of real systems, where parti- 
cles travel through non-repeating impurity environments. 
Also, the renormalization group predictions are predi- 
cated upon uncorrelated disorder, so that comparisons 
are not possible. 

Clearly, this predicament presents a serious drawback 
for a baseline study in novel territory. We are thus mo- 
tivated to undertake the first case study which is free 
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from this difficulty, by investigating the efFects of frozen 
disorder on a NESS without global current. In the follow- 
ing, we first describe our microscopie model and derive 
the associated Langevin equation. Next, we compute its 
universal scaling properties to two-loop order in an ex- 
pansion around the upper critical dimension d c = 5. To 
illustrate our predictions, we consider the structure fac- 
tor, since the latter is easily measured in simulations. We 
conclude with some comments and open questions. 

In the absence of quenched disorder, our model |ï(]] 
consists of particles and holes, distributed over a fully 
periodic regular lattice, with the usual nearest-neighbor 
attractive Ising interaction. The drive E = E(x, t)ê 
acts along a specific lattice direction ê, labelled 'paral- 
lel', and is characterized by an annealed random ampli- 
tude, E(x,t). Partiele- hole exchanges in the transverse 
directions occur according to the usual Ising energetics, 
while partiele hops along (against) E are enhanced (sup- 
pressed) by a factor of exp(+£/T) ( exp(-£/T)). For 
each parallel exchange, occurring at time t and lattice 
site x, a new value of E(x, t) is chosen from a symmetrie 
distribution with zero mean, so that no global current 
through the system is induced. Otherwise, the detailed 
form is irrelevant for universal critical behavior. 

Like its uniformly driven counterpart, this randomly 
driven system exhibits a continuous transition, from a 
disordered into an strip-like ordered phase. It is in 
the same universality class as an Ising lattice gas with 
Kawasaki dynamics, controlled by two temperatures 
Its scaling properties have been calculated to two-loop 
order in an expansion around its an upper critical di- 
mension d c — 3 jï(^[nj, and agree well with results from 
simulations |Ï2| . Intriguingly, the fixed point of this non- 
equilibrium system is Hamiltonian pij ], formally identi- 
cal to the one controlling the equilibrium scaling behavior 
of uniaxial ferromagnets with long-range dipolar interac- 
tions H. 

Next, we introducé frozen disorder into the dynamics 
of the system, by adding a quenched random contribu- 
tion, A(x), to the amplitude of the drive. Thus, we fo- 
cus on the non-equilibrium aspects of the problem, in 
particular since the statics of uniaxial ferromagnets with 
quenched disorder has already been discussed elsewhere 
Q . Writing Ê = [E(x,t) + A(x)]ê, we assume that the 
quenched component A(x) is distributed symmetrically 
around zero mean, so that no global current is induced, 
with second moment (A(x)A(x')) = 2au S(x — x'). Note 
that there is no (5-function in time here, consistent with 
the quenched character of this term. Again, the detailed 
form of the distribution is irrelevant for leading singular 
behavior. This type of randomness is easily implemented 
in a Monte Carlo simulation, where it corresponds to a 
time-independent local bias resident on each longitudi- 
nal bond. In the absence of A, the system reduces to the 
randomly driven (or two-temperature) lattice gas which 
we will refer to as the 'pure' system. 



The only slow variable of our theory is the local magne- 
tization, (f>(x, t). Since it is conserved, the effective equa- 
tion of motion for its long-wavelength, long-time behav- 
ior takes the form of a continuity equation, d t <p — — V ■ j. 
The current j here is composed of two terms, a deter- 
ministic and a noisy part. In the absence of E, it is 
given by Model B of the Halpcrin-Hohenberg scheme 
Hl, i.e., j = -V{5H/84>) + (, where H[<j>] is the usual 
Landau-Ginzburg-Wilson Hamiltonian. ( represents the 
thermal noise, modeling the relevant effects of the fast 
microscopie degrees of freedom after coarse-graining. It 
is Gaussian distributed with zero mean and second mo- 
ment ((i(x,t)(j(x,t')) = 2XNij5(x - x')S(t - t'). Here, 
i,j label the cartesian coordinates of C, and Nij is a 
noise correlation matrix. In equilibrium, the fluctuation- 
dissipation theorem (FDT) (!(]] requires iVy oc A 
non-vanishing E has two major effects on the equation 
of motion: First, it induces an additional contribution, 
je, to the deterministic part of j. Familiar from the uni- 
formly driven system, je = c((p)E is an Ohmic term, 
with the conductivity c(4>) cx 1 — <fi 2 reflecting the van- 
ishing of je in a completely filled or empty system. Sec- 
ond, since E singles out a preferred direction, j will bc- 
come anisotropic, so that all couplings have to be re- 
placed by appropriate anisotropic ones |ÏÖ1] . In particu- 
lar, there will be two diffusion coefficients, th and r±. 
Given that E acts like an extra noise along the paral- 
lel direction, th is enhanced over tj_, so that criticality 
is associated with tj_ — > while T|| remains positivo. 
Moreover, while still diagonal, the noise correlation ma- 
trix is no longer proportional to the unit matrix. Instead, 
C(x,t) splits into a parallel component, G , and a (d— 1)- 
dimensional transverse component, distributed ac- 

so that 
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the FDT is broken. Finally, we return to the divergence 
of the Ohmic current, V- je oc V- (1 — (f> 2 )E , and insert 

E = [E(x,t) + A(x)]ê. Letting V (d) denote the gradiënt 
in the transverse (longitudinal) directions, we recognize 
that the term dE(x, t) simply enhances the thermal noise 
dQ\ ■ In contrast, <9A(x) gives rise to a new term, cor- 
responding to a quenched random Langevin noise, which 
will be responsible for profound changes in the critical 
behavior, compared to the pure system. The remaining 
term, dcj) 2 (E + A), generates only couplings that are ir- 
relevant in the renormalization group sense. 

The resulting Langevin equation of motion must, in 
principle, be solved for <p(x,t), so that thermodynamic 
observables, such as correlation functions, can be com- 
puted by averaging over the thermal noise (, the annealed 
drive E(x, t) and the frozen disorder A(x). In practice, it 
is easier to introducé a Martin-Siggia Rose JÏ7| response 
field 4>, and recast the Langevin equation as a dynamic 
functional |19| , J[4>, 4>]. In this form, the averages over 
£, E(x, t) and A(x) are easily performed, and correla- 
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tion and response functions follow as functional averages weighted by exp(— 3)- Neglecting irrelevant terms, 
3 [<j>, $\ = ( d d x j dt4> {<j> - X [(rj_ - V 2 )V 2 + t,|Ö 2 + |v 2 3 ] } 



+ / d d x{ / dt\4>(x,t){n\\d 2 + n ± V 2 )4>(x,t) + dt dt'X 2 ^{x, t)(a {l d 2 - <tiV 4 )^(i, /ƒ)} 



(1) 



The term in the square brackets is easily recognized 
as the deterministic part of V ■ j. The two contributions 
that are quadratic in <fi arise from the averages over the 
thermal noise C,{x, t) and the quenched disorder A(x), re- 
spectively. Since £ is J-correlated in both space and time, 
the resulting operator is fully local. In contrast, the oper- 
ator associated with the quenched disorder is delocalized 
in time, reflecting the infinite correlation time of A(x). 
We emphasize that it consists of two terms: In addition 
to the coupling <7ii , arising from the average over A(cc), a 
second time-delocalized coupling <j± is generated under 
the renormalization group. It is rcquired for closure and 
must be included in (1), even though it did not appear in 
our original Langevin equation. Finally, we note that all 
coefficients here are functions of the microscopie model 
parameters, but - in the spirit of Landau theory - their 
explicit dependence need not be known. 

The first step towards the renormalization group anal- 
ysis of our theory consists in dimensional analysis. In- 
troducing a typical scale /j, to measure transverse mo- 
menta k±, we first note that, near criticality, parallel mo- 
menta scale as k\\ oc fï 2 . Next, requiring that all terms 
in (1) be dimensionless, we find Xt cx fi , so that the 
thermal noise becomes irrelevant and may be neglected, 
compared to the noise term associated with the frozen 
disorder. Furthermore, choosing the scale of (f> such that 
a±_ oc we find </> oc /j( d ~ 3 )/ 2 and (f> oc /j( d+5 )/ 2 , imply- 
ing that g oc /i 5 ~ d . The upper critical dimension follows 
as d c — 5, which is shifted upwards from pure system 
where d c = 4. Thus, our first major result emerges, 
namely, that the disordered system belongs to a different 
universality class than the pure system. Similar dimen- 
sional shifts have been observed in the uniformly driven 
system when frozen disorder is included [frUs!. However, 
we note that the e-expansions of the disordered and the 
pure systems do not coincide here, in contrast to random- 
field ferromagnets in equilibrium |Ï8| . 

Following Standard methods |^Ö| , we now consider the 
critical theory with insertions of Xr±(j)V 2 (f>, to two-loop 
order in e = d — d c . The bare correlation and re- 
sponse propagators take the anisotropic forms dictated 
by (1), which differ from their counterparts in Model B 
or the pure system. In particular, the correlator S(k,uj) 
becomes proportional to <5(w), by virtue of the time- 
delocalized noise. Noting that (1) is invariant under the 
anisotropic scale transformation ieii — > ax\\, xj_ — > xj_, 
4> — > <j> — > (a/3) -1 </), the effective dimensionless ex- 



pansion parameters of the theory can be identified as 
V = a±_T^ 1/2 g, and w = <r\\/ (a±_T\\). 

Next, we proceed to compute the one-particle irre- 
ducible vertex functions with n (n) external (f>- (</>-) legs 
and m insertions, denoted by in a doublé expan- 

sion in v and w. As in the pure system, r^} , T^, as 
well as are ultraviolet divergent and can be ren- 
dered finite by renormalizing (j), g, and tj_, respectively. 
In addition, however, r^g o-lso requires renormalization 
here, by redefining <r±. Beyond these four, no other in- 
dependent renormalizations are needed. For e > 0, we 
find exactly one stable, nontrivial fixed point (v*,w*), 
which controls the critical behavior of the theory below 
five dimensions. Surprisingly, the coupling w is found 
to be irrelevant under the RG, approaching its fixed 
point value w* — quite slowly, with a correction-to- 
scaling exponent lj w — ^e 2 {ln | — i|} + 0(e 3 ). v, 
on the other hand, flows much faster to a finite value 
v* = |e + (|) 3 Ce 2 + 0(e 3 ), controlled by a correction- 
to-scaling exponent u> v = e — (|) 2 Ce 2 + 0(e 3 ). Here, 
C = pp + | ln 2 — || ln 3 is a numerical constant. 

Incorporating dimensional analysis and the scale in- 
variance of the theory, the solution of the renormalization 
group equations at the fixed point yields the full scaling 
behavior of, e.g., the steady-state structure factor, 

S(k,t) = M - 4 +"S(A;||/ M 1 + A ,fe ± / M ,^;T±/M 1/l ') 

which serves to define the four indices v, rj, z, 
A f 



(2) 
and 



|l|. In particular, we find that v = ' 
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1 -- 7 i -31 ln | + 2 ln 4] + 0(e 3 ) is determined by the 
renormalization of t±, resulting in an exponent that is 
distinct from its counterpart in the pure system |ïïf| . The 
remaining two independent renormalizations, of the field 
<f> and the coupling <j± , give rise to two independent expo- 
nents, T) = -^f£ 2 + 0(e 3 ) and a = -^s 2 [11 + 31n|] + 
0(e 3 ). Scaling laws express the three exponents 77, z, 
and A in terms of these two quantitics. Specifically, 
the dynamic exponent z follows as z = 4 — 77 while the 
strong anisotropy exponent A is given by A = 1 — 77/2. 
Thus, we recover a scaling law, z = 4 — 2(1 — A), which 
is also obeyed by the pure system. In contrast, how- 
ever, the exponent 77 is not simply related to fj: Unlike 
the pure system, where 77 = 77, we find 7/ = 2 — ër = 
|j [| — 41n|] + 0(e 3 ) by virtue of the additional inde- 
pendent renormalization of the coupling a± . 

Other scaling laws determine the remaining anisotropic 
exponents pj. Some of these involve only A, and hence 77, 
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such as the two different correlation length indices v± = v 
and v\\ = v{\ + A), or the two different dynamic expo- 
nents, z± = z and z || = z/(l + A). Since A > 0, we 
find z\\ < z± so that transverse critical fluctuations de- 
cay more slowly than parallel ones. On the other hand, 
the four 77-like exponents, characterizing the anisotropic 
power-law behavior of the structure factor in real (RS) 
and momentum (MS) space, are determined by both r\ 
andfj: nf s = i], rjf s = (r) + 2A)/(1 + A), rjf s = ij + A, 
and r]^ s = [rj—{d — 3) A]/(l + A). Finally, we quote the 
order parameter exponent (3 = hv(d — 4 + 77 + A) which 
follows from the scaling form of the equation of state [^Ö) . 

In summary, we have analyzed the critical behavior of 
a randomly driven stochastic lattice gas, under the ad- 
ditional influence of a frozen local bias, acting along the 
parallel direction. This bias produces a time-delocalized 
Langevin current, with components in both, parallel and 
transverse, directions, that dominates the usual thermal 
noise. As a result, the critical dimension is shifted up- 
wards to d c = 5, and the critical behavior falls into a new 
universality class, distinct from the pure system. 

We conclude with two remarks. First, we return briefly 
to the microscopie mechanisms giving rise to the coupling 
ct_l, i.e., the transverse part of the frozen Langevin cur- 
rent. As discussed above, it is generated by a quenched 
bias acting along parallel bonds. Intriguingly, however, 
the same term also arises in a very different microscopie 
context, namely, if the quenched disorder resides in the 
nearest-neighbor hopping rates ||, modelled, e.g., by 
a Gaussian distribution of local barriers between sites. 
Characterizing this landscape by a coarse-grained poten- 
tial 0] U(x), the Langevin current takes the form of a 
potential gradiënt, VU(x). Averaging over U then leads 
to ƒ dt ƒ dt'X 2 (i>(x,t)(a±V 4 + K\\d 4 )4>{x 1 t') in (1) instead, 
with the latter term being irrelevant for critical behavior 
near d c = 5. Thus, two rather different microscopie mod- 
els, distinguished by the implementation of the frozen 
disorder, are described by the same mesoscopic Langevin 
equation, so that they both fall into the same univer- 
sality class. On the other hand, a quenched bias acting 
along transverse bonds gives rise to different universal be- 
havior |2Ï| . Work is in progress to clarify the intriguing 
connections between different microscopics with frozen 
disorder and their surprising degree of universality. Sec- 
ond, we wish to comment on the observability of our 
predictions in a computer experiment. To measure the 
full time-dependent structure factor, the leading singular 
part, captured by (2), must be separated from correc- 
tions to scaling, generated by irrelevant couplings such 
as, e.g., the thermal noise terms. Since the latter decay 
exponentially in time, only the former survives for large 
t. lts prime signature is its strong fc^ 4+ï? divergence as 
k± — > 0. Clearly, it would be very interesting to test our 
predictions with simulations. 
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